Higher Harmonics in Non-Linear Vacuum from QED Effects 
Without Low Mass Intermediate Particles 
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We show that in the presence of a slowly rotating strong transverse magnetic field there is an 
infinite spectrum of harmonic wave functions A„ due to the first order QED correction (in a^) 
given by the Euler-Heisenberg Lagrangian. The frequency shifts are integer multiples ±ul>o n of 
the magnetic field angular frequency rotation ujo = 2-nVrn and the several modes n are coupled to 
the nearest harmonics nil. This is a new effect due to QED vacuum fiuctuations, not exploit 
before, that can explain, both qualitatively and quantitatively, the recent experimental results of 
the PVLAS collaboration without the need of a low mass intermediate particle, hence may dismiss 
the recent claim of the discovery of the axion. 



In quantum field theory the vacuum is not empty, in 
the sense that virtual particle pair creation and annihi- 
lation are taking place. Considering the case of QED 
-Quantum Electrodynamics- in the presence of external 
electromagnetic fields the vacuum shows properties close 
to the properties of any other optical medium, as dichro- 
ism and birefringence. These effects can be considered as 
first order corrections to the Maxwell theory [l], 0, |3, 01 , 
and the basic diagrams, correcting the photon propaga- 
tor, are shown in figure [1] In the presence of low mass 
neutral particles coupling to two photons, like the tt^ and, 
more interesting, the axion -Primakov effect- or the gravi- 
ton, there are additional non- linear effects (see figureTT]), 
which may be useful in directly detect the axion jg, l6| . 

In this letter we present and exploit an effect not con- 
sidered before, in vacuum a strong rotating magnetic field 
generates an infinite set of higher harmonics. We show 
that recent experimental results 01 can be explained in 
the basis of pure QED which dismisses the recent claim 
of the discovery of the axion . 

In the QED lowest order correction (radiative cor- 
rections in a^) to the Maxwell Lagrangian £Maxwoii = 
—Ffj^^F^'^ is given by the regularized Euler-Heisenberg 
Lagrangian 



£(2) ^ ^ U iF^.F^-f + 7 (e^^'PF^^Fsp) 



(1) 



where ^ = 2a/(45i?^) ^ o? and the critical field is 
Be = mlc^lieh) = 4.414 x 10^ T Q. The usual pro- 
cedure is to take a decomposition of the gauge connec- 
tion F into an external component Fq plus an internal 
component /^^ = dy^A^ — d^A^, such that we obtain 
Fyv = i^o,/i!^+/Ai!^j a-iid consider the lowest order diagrams 
of photon interactions due to the background electromag- 
netic fields. There are two main processes that may be 
considered, 'photon splitting and photon dispersion. De- 
pending on the practical application and experimental 
setups there are other processes, such as, photon-photon 
interactions fsl and other photon alternative processes [oj 
that are relevant. In this work we address only the pro- 
cess of photon dispersion in the low frequency and weak 
field limit, i.e. oj <C 2meC^/h and |F| ^ mec'^/{he). 
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FIG. 1: Feynman diagrams corresponding to the electron loop 
that contributes to the photon dispersion effect and the axion 
exchange. 



The relevant diagrams are the ones containing one inter- 
nal fermion loop with two external photons and an even 
number of exchanged photons between the fermion loop 
and the external electromagnetic field (see figure [l}. 

We obtain that the relevant radiative corrections to 
the usual classical wave equation in order is linear in 
the photon field A [i| 



(V^ - \d^)A - A.A + 0{a*) , 
where the matrix A has the eigenvalues [1, 0| 



(2) 



A|| = 14^(t^/c)2|Q|2 and = 8(,{uj/c)'^\Q\^, (3) 

where Q = k x Eo/|k| + k x (k x Bo)/|kp and k stands 
for the light wave vector. Eo and Bo are the external 
fields and the directions || and ± correspond respectively 
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to the parallel and transverse directions to the external 
magnetic field. 

We consider linear polarized light traveling in vacuum 
in the z direction with wave number k and frequency uj 
under a slowly rotating transverse magnetic field rotating 
with angular frequency cjq = 'i.iiVm <C oj such that we 
have A = [A^.Ay.Q), k = (0,0, /c), w = ck, Bo{t) = 
Bo[sm{cjot),cos{uJot),0], Eq = and |Q| = cBq. In 
these conditions the gauge field depends only on z and we 
can write the equation ^ for A in terms of the parallel 
direction (||) and transverse direction (_L) as derived by 
Adler However these directions are rotating with an 
angular frequency loq, therefore we take the directions of 



the decomposition of the fields to be the directions of 
the magnetic field at some initial time t = 0. Then the 
eigenvalues of the matrix A become time dependent since 
they are valid with respect to the magnetic field direction 
such that we have 



A|| (t) = A|| cos(ijJo t) — X± sin(u;o t) 
A±{t) — A|| sin(a;o t) + A± cos(a;o t) 



(4) 



with A|| and A_l given in ([3]). Considering the above eigen- 
values we obtain the field equation 



J 



id'. 



; AO ) - (A|| cos(c^n t) - \i_ sin^ t)) (Ajj, 0) + (A,, sin(cc>o t) + A^ cos(coo t)) (0, 



(5) 



r 



where the Ajj and correspond to the decomposition of 
the field along the directions parallel (||) and transverse 
(_L) to the magnetic field at the initial time t — 0. For 
constant magnetic field {luq = 0) these equations imply 
that the vacuum acquires a birefringence due to the ex- 
istence of two different refractive indices = I — Xi/2 
along both the directions i =||,-L. For slowly rotating 
magnetic fields the usual procedure is to take the zeroth 
order expansion of the cosines and sines that hold the 
same effect [6] . However we show next that even for small 
ojQ equation ([5]) implies, not just a vacuum birefringence, 
but instead a coupling between the neighbor harmonics 
separated by frequency shifts of iwo- It is important to 
note that equation ([5]) is valid for slow rotating external 
fields (see [4] for further details), i.e. luq ^ cu. 

We can cast equation ([5]) in complex form by using the 
definitions X± — A|| ±iA_L, A± = A|| ±1^4^ and the usual 
relations for the trigonometric functions and consider a 
mode decomposition of the form 



A±{z,t)=Y: 



c fc„ = w„ = cj + t^o , = { ±-j n , (6) 



Xf^+ X'^j_ , Ox = arctan ^ ^ 



Given the above mode decomposition the differential 
equation becomes a recursion relation in the several fre- 
quency modes An 



V2cXa 



An — An+1 + An-1 



(7) 



where An = = e ^'"^^An .We neglected a second 
derivative term An and considered a dispersion relation 



that renders a tower of refractive indices in each direc- 
tion. This is the same approach of [4] in consistent field 
approach. Here we have for each mode n and each direc- 
tion i =11, -L a different refractive index 



Wo 



kc = ii'l{Lu)u; , n^iLu) = 1 — n . 



(8) 



The solutions up to an overall phase for the recursive 
relation ([7|) are modified Bessel function of the first kind 



A„(<) = (-!)"/„ (zr) 



V2cXot 



(9) 



We note that indeed the second derivative of An is of 
order a* and can be safely neglected. These functions are 
normalized for any r and for even n the functions are pure 
reals while for odd n the functions are pure imaginary. 
The only mode that is non null at t = is loiir) such 
that only for larger values of t the amplitudes of the other 
modes increase. These properties are important when 
imposing the initial conditions, i.e. matching our solution 
with the incident wave at the initial time t = 0. 

The results obtain so far are very interesting and 
describe both vacuum birefringence and a generalized 
dichroism in the presence of a strong rotating magnetic 
field. The tower of refractive indices ^ correspond to 
a birefringence effect quantized in terms of the quantity 
ojq/uj. The Bessel functions solutions account for a gen- 
eralized dichroism effect that decreases the amplitude of 
the incident wave in one mode and increases it in the 
nearest modes, we note that the QED effects considered 
here do preserve momenta such that there is no real loss 
of energy, it simply is transfered from one mode to the 
nearest modes. These effects can be understand as fol- 
lows, from equation ([7]) we see that the variation of each 
mode amplitude An over time contributes to the nearest 
modes A„±i, this accounts for the generalized dichro- 
ism effect in the sense that no real absorption happens. 
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therefore we have a true rotation of polarization (in the 
sense considered in 0). As for birefringence we note 
that by combining the nearest modes (e.g. n|| = n and 
n± — nil) in different directions we will have a net 
birefringence effect given by n_L ~ n|| — ±ujq/uj such that 
each mode combination correspond to a wave with elliptic 
polarization. It is also interesting to note that, although 
the description for the propagating wave is given in terms 



of the collective field A, we may as well interpret these 
solutions as cross section for photons. In such interpre- 
tation the allowed configurations for a single photon are 
independent in both spatial dimensions (||,-L), such that 
indeed we have elliptic polarized waves. 

A solution for the electric field compatible with lin- 
earized polarized waves at the initial time r = is 



+ 00 



^ = Eo\ {hn{iT)coa{kz-Lj2nt + O2n)-ihn-i{iT)sm{kz-uj2n-it-O2n-i)]j{sin{0o),cos{eo)). (10) 

I 



L n— — oo 



Here Eq is the electric field amplitude, ujn — u+uo n and 
6n = —{S\ + tt/^) n. This field is defined in the referential 
of the parallel and transverse direction to the magnetic 
field at the initial time t = such that 9q stands for the 
angle between the original linear polarization direction 
and the magnetic field direction at t = 0. 

Take the n|| — ni mode in one of the directions to- 
gether with one of the other modes n± = n2 > ni in 
the other direction such that this combination has now 
an elliptic polarization and the difference of argument 
of the trigonometric functions corresponding to each di- 
rection correspond to a phase difference of A(p„2_„j = 
(n2-ni)(wot + 6'A + 7r/4). 

We proceed to compare our results with the recent ex- 
perimental results of the PVLAS collaboration \f\. For 
the apparatus of that experiment an incident polarized 
laser beam traverses a region under a strong rotating 
magnetic field such that it is possible to measure po- 
larization rotations (that correspond to vacuum dichro- 
ism) and ellipticities (that correspond to vacuum bire- 
fringence) . 

We can write the amplitude measured in the experience 
for our solutions as 



+ CKj 



2 J2 \In{iT)\ cos(n (wq t + + ^0 + ^/4)) 



+r]o cos(ws t + 9s)+ To 



(11) 

Here 770 cos(ci;5 1 + 9s) stands for the modulator wave 
and heterodyne detection is used to identify the several 
harmonics present due to rotations. r„ are corrections 
that encode uncompensated dispersions for the several 
modes n due, either to systematic deviations, or unac- 
counted processes. In (|lip we sum over rotations in one 
direction only, that is why we have a factor of 2 in front of 
the sum. By expanding the square of this expression we 
obtain the angular frequency spectrum and the respec- 
tive amplitudes and phases for each of them. We list the 
more significant frequencies that can be measured near 
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es±n (Ox + 60+ 7r/4) 



TABLE I: The relevant angular frequency spectrum measured 
near the modulator frequency lus with respective amplitudes 
and phases. There is an infinite set of other frequencies that 
have lower amplitudes. 



the modulator frequency 0^5 and list them in table HI 

We take the values of figure 2 of the original refer- 
ence [tI (expressed in I.S. units) 



w = 1.772 X 10+15 s"i 



Az = cAt = 4.6 X 10+^^ m 



LUq 



(12) 



Bn = 5.5 T 



The theoretical results obtained above are valid for these 
experimental conditions, in particular the magnetic field 
is slowly rotating (loq <C w)and obey the low frequency 
and weak field approximations, respectively uj ~ 10^^ ^ 
2me c^/h - 10^1 and - 10 < c^/{eh) ~ lO^^. For 
the values given in (|12p and using the finesse of the Fabry- 
Perot cavity « 7 x 10^ (corresponding to an increase 
of optical path N « 4.46 x 10"*) ^ we have from ^ 
with r = 0.0088 :'3:] that |/o(ir)| « 0.9998, |/i(ir)| « 
4.384 X 10-^ |/2(jT)| w 9.609 x lO'^ and Ihiir)] « 
1.404 X 10"^. We use these values to fit our solutions 
to the data expressed in figure 2 of reference Q. We 
considered that the background noise is « —135 dBVrms 
and assumed values for rjQ and the F^'s that allow us to 
fit the experimental results. These values are given as 
an example and are not based in a detailed analysis. In 
table ini we list a possible set of values (in dBVrms) that 
fits the experimental results of Q . We note that the value 
obtained for the mode n = 3 (« —145.65 dBVrms) is 
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frequency 
[Hz) 


T 

-i-cxp 

{dBVrms) 


Voln 

{dBVrms) 




506 


^ —51 




~ +16.93 


506 ± 0.35 


^ -91 


« -90.70 


+0.704 


506 ±0.70 


f=s -118 


-117.30 


^ -0.704 


506 ± 1.05 


< -135 


^ -145.65 





TABLE II: We list the experimental values taken from figure 2 
of reference [3] with the theoretical values given by the /„(ir) 
up to n = 3. We considered 770 « —68.63 dBVrms. This 
value, as well as the value of the r„'s, are only given as an 
example that it is possible to fit the data to our solutions. 
The theoretical value found for the amplitude of a rotation 
corresponding to n = 3 is ~ —145.65 dBVrms and is therefore 
below the noise level of the experiment that we considered 
(« -135 dBVrms). 

below the noise level (« —135 dBVrms), hence it could 
not be detected for the available data. 

To conclude we compute the expected rotation 
values corresponding to the second harmonic n = 2. 
By direct geometrical analysis we conclude that 
the angular rotations of polarization are given by 

the generic expression ip — 9q ~ arctan ^ ^^ut ^ , 

where E°"* = {E^^'^\ Ef^) stands for the radiation 
electric field reaching the detector. For a given har- 
monic n we have the possible combinations E°"*+ = 
[/o sin(0o) cos((5(f, z)), /„ cos(6'o) cos((5(t, z) + ojq At + 
arctan(7/4) + 7r/4)] and E™'" = [/„ cos(6'o) cos{5{t, z) + 
Wo At + arctan(7/4) + 7r/4), /q sin(6lo) cos((5(t, z))]. Here 
the phase 5{t,z) = kz — tot is the phase of the wave 
when reaching the detector and the sine and cosine of 
^0 is due to the angle between the original polarization 
and the magnetic field when the radiation enters the 
magnetic field as expressed in the solution pop . 

Based in the values p2)l with 6q = tt/A (that corre- 



sponds to the maximum effect) and assuming S{t, z) = 
we obtain for the rotations corresponding to n—2 f^^2 ~ 
3.78 X 10^^ rad. This value is one order of magnitude 
higher than the experimental values measured (~ 10"''). 
A possible explanation for this result may be due to a dif- 
ferent initial value of the angle between the polarization 
and the magnetic field, ^o- If we consider for instance 
9q — 7r/48 we obtain iPn=2 ~ 2.48 x 10~^ rac? which is 
closer to the experimental observed rotation. Other pos- 
sibilities may be due to an effective lower path length due 
either to decoherence or unaccounted losses. 

So we have presented a new effect of vacuum polariza- 
tion rotation due to external strong magnetic fields based 
in pure QED. This effect generates an infinite number of 
harmonics and explains, both qualitatively and quantita- 
tively, the PVLAS experimental results . These exper- 
imental results show that the PVLAS experiment gives 
(as far as we know for the first time) a beautiful demon- 
stration of the QED properties of vacuum and have be- 
fore hand registered an effect that have not been justified 
at theoretical level. Our results may dismiss the recent 
claim of the finding of the axion [l^ . Both for future use 
in experiments currently being planned and for a correct 
interpretation of already detected (and to be detected) 
physical phenomena [ll[, a full analysis of existent data 
is required. In particular, we note that, if the noise level 
can be decreased, the third peak {n — ±3) should be 
detected. 

Erratum 

A new mechanism was proposed to explain the ob- 
served spectrum in the PVLAS experiment 0]. We have 
predicted a series of new signals with frequencies uj-\-nujQ, 
for n integer, where uj is the laser frequency and the 
angular frequency of the rotating magnetic field. But two 
mistakes were made in this work. First, the secondary 
signals have frequencies uj + 2noJo. Second, our numerical 
estimates are wrong by a factor of 10^ [l^, and therefore 
cannot explain the observed signal. It should be stressed 
that signals at frequencies uj + nujo can still be generated, 
but they are not relevant to the PVLAS experiment. 
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